Spin fluctuations and pseudogap in the two-dimensional half-filled Hubbard model at 

weak coupling 
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Starting from the Hubbard model in the weak-coupling limit, we derive a spin-fermion model 
where the collective spin excitations are described by a non-linear sigma model. This result is 
used to compute the fermion spectral function y4(k,tj) in the low-temperature regime where the 
antiferromagnetic (AF) coherence length is exponentially large ( "renormalized classical" regime). 
At the Fermi level, A(\^F,(JJ) exhibits two peaks around ±Ao (with Ao the mean-field gap), which 
are precursors of the zero-temperature AF bands, separated by a pseudogap. 

PACS Numbers: 71.10.Fd, 75.10.Lp 



X3 
O 

(N 
> 
oo 
m 



o 



o 



X 
J3 



I. INTRODUCTION 



In the last two decades, the discovery of heavy-fermion 
compounds, high- Tc superconductors and organic con- 
ductors has revived interest in strongly correlated elec- 
tron systems. Of particular interest are metallic phases 
which, although conducting, are not described by Lan- 
dau's Fermi liquid theory because of the absence of well- 
defined quasi-particle excitations. A well-known example 
is given by the normal phase of high- Tc superconductors. 
Instead of quasi-particles, these systems exhibit a pseru 
dogap at low energy as shown by many experiments. □ 
Although the origin of the pseudogap is still under de- 
bate, it is generally believed that antiferromagnetic (AF) 
fluctuations play a crucial role. 

In this paper, we consider the pseudogap issue on the 
basis of the half-filled 2D Hubbard model. We consider 
only the weak-coupling limit U <C i (?7 is the local 
Coulomb repulsion and t the intersite hopping ampli- 
tude). [In the strong-coupling limit at half-filling, the fi- 
nite temperature paramagnetic phase is a Mott-Hubbard 
insulator with a (charge) gap of order U . At T = 
there is a transition to a Neel antiferromagnetic state 
Although the ground-state is AF, long-range order is de- 
stroyed by classical fiuctuations at any finite tempera- 
ture, in agreement with the Mermin- Wagner theorem. 
Nevertheless, below a crossover temperature Tx (of the 
order of the mean-field transition temperature) , the sys- 
tem enters a renormalized classical regime where AF cor- 
relations start to grow exponentially. Contrary to the 
3D case, at the zero-temperature 2D phase transition 
the system goes directly into the (Neel) ordered state 
where the fermion spectral function A(k, uj) exhibits two 
well-defined quasi-particle (QP) peaks corresponding to 
the Bogoliubov QP's. By continuity, the two-peak struc- 
ture in A(k, oj) cannot disappear as soon as we raise the 
temperature. As pointed out in Ref. ^, the only possible 
scenario is that at finite but low temperature the fermion 
spectral function exhibits two (broadened) peaks which 
are precursors of the T = Bogoliubov QP's, separated 



by a pseudogap. We therefore expect the presence of a 
pseudogap at finite temperature, due to the strong (clas- 
sical) AF fluctuations. 

Clearly, traditional mean-field techniques fail to de- 
scribe these phenomena. For instance, the random-phase 
approximation (RPA) predicts a finite temperature phase 
transition which is forbidden in 2D by the Mermin- 
Wagner theorem. More sophisticated approaches are 
therefore required. In the weak-coupling limit, the pseu- 
dogap formation has been considered within the fiuc- 
tuation exchange (FLEX) approximation: .and the two- 
particle-self-consistent (TPSC) theoryoO'B which both 
satisfy the Mermin- Wagner theorem. Only the TPSC 
theory predicts the formation of a pseudogap in the 
fermion spectral function A(k, lo) at low temperature. 

The aim of this paper is to describe an alternative ap- 
proach to the 2D half-filled Hubbard model in the weak- 
coupling limit. We first derive a spin-fermion model 
where the collective spin excitations are described by a 
non-linear sigma model (NLctM). The spin- wave velocity 
and the coupling constant of the NLuM are expressed in 
terms of the ground-state properties of the system. Solv- 
ing the NLctM in a "large- A/"" limit, we then compute the 
fermion spectral function A(k, uo) to lowest order in the 
spin-fermion interaction. At the Fermi level, A(ki?,w) 
exhibits two peaks around ±Ao (with Aq the mean-field 
gap) which are precursors of the zero-temperature AF 
bands, separated by a pseudogap. We compare our re- 
sults with those of the TPSC theory. 



II. MODEL 

The two-dimensional Hubbard model is defined by the 
Hamiltonian 

H = -t (ct.Cr'. + h.C.) + {/ ^ rirT^r^ (2.1) 

{r,r'),(T r 

where t is the intersite hopping amplitude and U the on- 
site Coulomb repulsion. Cra is a fermionic operator for 
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a cr-spin particle at site r [a =t,i), and n-ca = cl^Cra- 
(r, r') denotes nearest-neighbor sites. We take the lattice 
spacing equal to unity and h = ks = I throughout the 
paper. 

Since spin fluctuations play a crucial role in the Hub- 
bard model at half-filling, it is convenient to intro- 
duce auxiliary fields describing these collective excita- 
tions. The standard approach is to write the interac- 
tion part of the Hamiltonian in terms of charge and 
spin fluctuations, i.e. rirfnrx = [(cjcr)'^ — (cJcrzCr)^]/4, 
and then perform a Hubbard-Stratonovich transforma- 
tion by means of two (real) auxiliary fields Ac and As 
[cr — (cr]-, Crx)"^]. Although this procedure recovers the 
standard mean-field (or Hartree-Fock) theory of the Neel 
state within a saddle-point approximation, it leads to 
a loss of spin-rotation invariance and does not allow 
to obtain the spin-wave excitations. Alternatively, one 
could write nriUri in a spin- rotation invariant form, e.g. 
iir^Tiri = — (cJ<TCr)^/6 whcrc (T deuotcs the Pauli matri- 
ces, and use a vector Hubbard-Stratonovich field. Such 
decompositions, however, do not rearoduce the mean- 
field results at the gaddle-point levehQ 

As noted earlier ,Q'la this difficulty can be circumvented 
by writing nrt^ri = [(cj^r)^ — (cJ<t • J7rCr)^]/4 where fir 
is an arbitrary unit vector. Spin-rotation invariance is 
maintained by averaging the partition function over all 
directions of Jlr- In a path integral formalism, Or be- 
comes a time-dependent variable. After the Hubbard- 
Stratonovich transformation, the partition function is 
given by Z = / ^[c^c]/ V[Ac, As,ft]e~^ with the ac- 
tion 

S = So + Y: /'dr{i(AL + AL) 



A. 



(2.2) 



So is the action in the absence of interaction. Since 
charge fluctuations are not critical (even when T ^ 0), 
they can be treated at the saddle point (i.e. Hartree- 
Fock) level. Their effect is to re norm alize the chemical 
potential fi from U/2 to 0. Eq. ( |2.2| ) then corresponds 
to a spin-fermion model where the fermions interact with 
their collective spin degrees of freedom (ArOr). [We now 
denote A^r by A^.] Below the crossover temperature Tx, 
i.e. when T <C Tx, low-energy excitations correspond 
to orientational spin fluctuations described by the unit 
vector field fir- We can then consider Ar within a sad- 
dle point approximation, i.e. Ar = Ao(— 1)"", where the 
fluctuations of Aq are ignored. In order to compute the 
fermion spectral function A(k, oj), one should flrst deter- 
mine the effective action S[ft] of the unit vector field ft. 



III. SPIN FLUCTUATIONS 

The effective action S[fl] is obtained by expanding 
around the Neel state. We first introduce a new field 



(f) defined by (f)r = i?Jcr, where Rr is a SU(2)/U(1) ma- 
trix which rotates the spin-quantization axis from z to 
ftr {RrCTzRl = fir ■ cr) . In tcrms of this new field, the 
action becomes 

fP 

S = S-MF + dT(f)lRldrRr(t)r 
r Jo 

-t V / dT[(l)l{RlRr. - l)(l)r' + C.C], (3.1) 



where Smf = Sq + ErfdriAl/U ~ Ao(-l)"-0t(T,0r). 
Within a saddle-point approximation with fir = z (i?r = 
1), i.e. ignoring spin fluctuations, one recovers the mean- 
field action Smf of the Neel state. The value of the order 
parameter, Aq = {U /2){—lY {(j)laz(j)r) , is obtained by 
minimizing the free energy. In the weak-coupling limit, 

this gives Aq ~ te"^''\/^i 

Low-energy spin excitations correspond to fluctuations 
of the unit vector field fi_|aM)und its saddle-point value. 
The standard proceduraijij'Ll is then to assume at least 
local AF order and write fir = nr(l - L^y/^ + (-l)''Lr, 
where the (Neel) order parameter field rir is slowly vary- 
ing in space and time and Lr is a small canting field 
(|nr| = 1, Lr • rir = and |Lr| <C 1). Integrating out 
both and L yields the action of the NLcrMBy In the 
strong-coupling limit U ^ t, one recovers the action de- 
rived from the Heisenberg model. 

As we now verify explicitly, the small canting field Lr 
gives negligible contributions to the parameters of the 
NLcrM in the weak-coupling limit [/ <C i. If we identify 
fir with the slowly varying Neel field, fir ~ rir, the ef- 
fective action ^[n] is readily obtained. Integrating out 
the fermions in Eq. (3.1) and taking the continuum limit 
in space, one obtains to lowest-order in gradient (i.e. in 
drR and VrRjB 



S[n] 
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dVdr[xT(a.n)2 + p^](Vrn)2 



(3.2) 



where x'± is the uniform transverse spin susceptibility 
in the mean-field state and p° = -{{K)mf/2 -I- n']_)/4 
the spin stiffness. Here {K)mf is the mean- value of the 
kinetic energy and the correlation function of the 
transverse spin current (j^ or j^). Eq. (3_^) should be 
supplemented with a short-distance cutoff (in momentum 
space) A ~ S^q^ , since short-range AF order cannot be de- 
fined at length scales smaller that the coherence length 
^0 t/Ao. Using the mean-field action Smf, one obtains 
{N is the number of lattice sites) 
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t^Ag sin^ kr, 

k 



(3.3) 
(3.4) 



2 I Ag)^/^ is the Bogoliubov quasi- 
particle excitation energy in the mean-field state (ck = 



where = (e^. 
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-2t(cos kx+cos ky) is the di spersion of the free fermions) 



We can verify that Eqs. 



can be directly ob- 
in the weak-coupling 



tained from the results of Refs. 
limit {U ^ t). The value of the spin- wave velocity 
c ~ v^pfTxT ~ t{U /tY^^ also agrees with the weak- 
coupling limit of the RPA result.Ej The approximation 
Sir ~ Hr is therefore justified when [/ <C While it re- 
stricts the validity of our approach to the weak-coupling 
limit, it makes the computation of fermionic correlation 
functions considerably simpler, sin ce t he fermions couple 
directly to the Neel field [see Eq. (|2^)]. 

We solve the NLctM within a "large- A/"" approach by 
extending the number of components of the uni t vector 
rir from 3 to A/". When A/" — > oo, the actipa (|]^) can be 
solved exactly by a saddle-point method.liJ Fig. 1 shows 
the resulting crossover diagram as a function of the di- 
mensionless coupling constant g = Ag ^ AcAf/ of the 
NLcrM. In the weak-coupling limit of the Hubbard model 



(U ^t),g^ cAoM/{p°t) (X 



is exponentially 



small. This implies that the ground state has AF long- 
range order with very weak quantum fluctuations. This 
magnetic order persists in the strong-coupling regime 

{U ^ t) wheiM g ^ (jc — 4:TT (see Fig. 1) in agreement 
with conclusions based on the Heisenberg model (for a 
square lattice). At finite temperature, magnetic long- 
range order is suppressed as required by the Mermin- 
Wagner theorem. The dominant fiuctuations are classi- 
cal since the gap m in the spin excitation spectrum (see 
below) is much smaller than the temperature (this regiaie 
is known as "renormalized classical" in the literaturella) . 

Since we are primarily interested in the fcrmion spec- 
tral function A(k, uj) at finite temperature, we shall con- 
sider the action S[n] in this regime. In the large- A/" limit, 
it reads 



S[n] 



A^ 
25c 



^(c.2 + cV+m2)|n(q,^c..)|^ (3.5) 



where we have introduced the Fourier transformed field 
n(q, iajjy) {lOi, is a bosonic Matsubara frequency). The 
length of the vector rir is no longer fixed to unity. In the 
large-A/" solution, the constraint jrirj = 1 is imposed only 
on average (via the Lagrange multiplier to).Ej The mass 
771 of the spin fluctuation propagator (a = 1 • • -Af) 



gc/U 



,2„2 



vn 



is determined by the saddle-point equation 

1 



1 = E 



(3.6) 



(3.7) 



In the renormalized classical regime, we can neglect quan- 
tum fluctuations. This approximation is excellent in the 
weak-coupling regime {U <^ t) since quantum fluctua- 
tions are weak {g <^ (jc, see Fig. 1). From Eq. (3.7), 



we then obtain the AF coherence length ^ = c/777 ^ 
A-^ exp{2Tr p^JAfT). 

Note that we ex pec t also a term m^luj^l/ujsi in the de- 
nominator in Eq. ( |3.6| ) . This term comes from the damp-, 
ing of spin fluctuations by gapless fermion excitations .113 
It is missed in our approach since we expand around the 
zero-temperature AF state which has only gapped quasi- 
particle excitations. Fluctuations are classical when 
m <^ T and Wsf ^ T. Both conditions are satisfied 
in the renormalized classical regime iTj-jjC Tx) since 
Wsf ~ (critical slowing down)J3Qll3 



IV. SPECTRAL FUNCTION 

Knowing the effective action S\r\\ of the spin excita- 
tions [Eq. ( |3.5| )], we are now in a position to compute 
the spectral function ^(k , uj) — — 7r~^ImG'(k, oj) from the 
spin-fermion model (2_^). Here G{k,uj) denotes the re- 
tarded part of the fermionic Green's function. By inte- 
grating first the fermions and then the spin fluctuations, 
we can write the Green's function as 

G(r - r', T - r') - 1 1 1?[n]e-^HG(r, r; r', r'ln). 

(4.1) 

G'(r, r; r', r'ln) is the Green's function for a given con- 
figuration of n: G^^[n] = Gq ^ -I- Ao(—l)'' cr • rir, where 
Go is the Green's function of the free fermions. Since 
^[n] is G aussian in the large-A/" limit, the averaging in 
Eq. ( [4.l| ) is easily done. The result can be written as 
G~"'^(k,za;„) = Gq ""^ (k, 7a;„) — S](k, 7w„) (wn is a fermionic 
Matsubara frequency). 

We consider the lowest-order contribution to the self 
energy E (Fig. 2): 

S(k, iw„) = Aq— E ■^x(q, ii^u)Go{\<i - Q - q, iujn - iuj„) 



^0-77 2^-2—1 



cN q- + 2 ^^^^ _ ek-Q-q 



(4.2) 



where the last line has been obtained in the classical limit 
{lu^ = 0) and Q = (tt, tt). At low temperature when 
^ 00, the sum over q in Eq. (4.2) diverges in 2D due 
to the contribution of long wavelengths (q ~ 0). We can 
therefore expand — ek-q-q = Ck-q — Ck — Vk • q around 
q = (vk is the velocity of the free fermions). Let us 
first consider a particle at the Fermi level. One easily 
finds that the imaginary part of the retarded self-energy 
{iuJn ^ OJ + zO+) takes the form 



E"(k 



P% 



th 



(4.3) 



where ^th = |vk|/T is the De Broglie thermal wavelength. 
Since ^ grows exponentially below Tx, it quickly becomes 
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larger than ^th- As a result, liniT^o /■Cth = oo and 
I]"(ki?,a; = 0) diverges at low temperature in contra- 
diction with the Fermi-liquid theory hypothesis. Thus, 
the lowest-order perturbation result shows that quasi- 
particles are suppressed by spin fluctuations when T <C 
Tx- This phenomenon is accompanied by the formation 
of a pseudogap. For \uj + ek| ^ |vk|/Ci the-real and 
imaginary parts of the self-energy are given bylHI 



E'(k,w) 



ek 



S"(k,o.) 



47r/90|w + ek| 



(4.4) 



Note that the condition + ek| ^ |vk|/^ is satisfied for 
any value of uj except in an exp onentially small window 
around ui — — ek- From Eq. (14), we deduce the spectral 
function 
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t: (l^2 _ Elf + 72 ' 



Aik, oj) exhibits two peaks at ±i?k that are precursors 
of the AF bands that exist in the T = ordered state. 
The width of these peaks is given by 7/A0 ~ TAo/p° ~ 
rp^-iTxyftJij ^ Xhe precursors of the AF bands are sepa- 
rated by a pseudogap. In particular ^(kj^jCj) vanishes 
at w = 0. 

WhenT^O (7-^0), 



^(''•")-5('-^t)'>-^'> + 5('-t) 



^(t^ + £;k), 

(4.6) 



temperature gapjl^ 
proportional to T.t 



which is the spectral functi on o f the T = AF state. 
Thus, the simple self-energy (4.2) predicts that the pseu- 
dogap evolves smoothly into the gap of the ground-state 
when T ^ 0. It should be noted that neglecting quan- 
tum fluctuations is justified only at low energy |w| < T. 
In particular, the precise location of the peaks around 
±Ao should depend on quantum fluctuations since Aq ~ 
Tx > T. 

The spectral function yl(k, ^fjlnffi'q- ( |4.5| )] is similar to 
the result of the TPSC theory.M In the latter, the po- 
sition of the majajaa. in A{k.p^uj) scales with the zero- 
\ i3,and the width of these two peaks is 
jlIi3 These two features agree with our 
conclusions. This similarity is not surprising sin ce i n both 
approaches a paramagnon-like self-energy [Eq. (4.2)] with 
a similar spin susceptibility [Eq. (3^)] is used to obtain 
the spectral function. The main difference comes from 
the spin fluctuation propagator x- While x comes from 
the NLcrM (which is itself based on an expansion around 
the ordered AF state), it is obtained by considering the 
paramagnetic phase in the TPSC theory. As a result, the 
basic parameters entering the spectral function A(k, to) 
[Eq. ( |4.5| )], namely the T = order parameter Aq and 
the T = spin stiffness Pg, do not appear in the TPSC 
theory. Instead, A(k, uS) is expressed only in terms of the 
paramagnetic properties of the system. 



Two comments are in order here. The validity of 
does not include vertex correction, may 
The importance of these corrections is a 



Eq. (^, whi 
be questioned. 

long-standing problem which is still under debate. Ver- 
tex corrections are expected to play a crucial role when 
higher-order self-energy contributions are taken into ac- 
count. The FLEX approximation, which sums up contri- 
butions to all order without vertex correction, does not 
predict the formation of a pseudogap in A(k, uj) at low 
temperaturcQ (see Ref. || for a detailed discussion of the 
FLEX approximation) 



In the spin-fermion model deflned by Eqs. (2.2) and 
( |3.2[ ), there are only two (transverse) spin excitation 
modes, as expected when only orientational fluctuations 
are important (T <C Tx). Unfortunately, thi s pr operty is 
lost in the large-TV limit of the NLcrM [Eq. (|]|)], where 
both transverse and amplitude fluctuations are allowed. 
Following Ref. A(k, lS) can be obtained exactly when 
^ — > 00 by summing all the self-energy diagrams. The 
result. 



A(k,w) 



3« 



27rAg 



(c^^-4)i/2(c. + ek)exp(-^ 



o 2 
6 u 



A2 



x[0(t^- |ek|)-0(-c^- |ek| 



(4.7) 



sho ws t wo broad incoherent features, located around 
±-y/2/3Ao for Ek = , instead of the correct T = 
limit given by Eq. (4.6). The correct limit is obtained 
only when, amplitude fluctuations are frozen in the limit 
^ cxd.E3Ej We therefore conclude that our approach, 
which is based on the large-A/" solution of the NLcrM, 
must break down at very low temperature. The fact that 
the spectral function A(k, w) derived from the lowest- 
order self-energy contribut ion does reproduce the correct 
result when T ^ [Eq. ( |4.5| )] appears somewhat acci- 
dental. A correct treatment of the T — > limit must 
freeze the amplitude fluctuations of the Neel field n. 



V. CONCLUSION 

We have described a new approach to the pseudogap 
in the half-filled 2D Hubbard model at weak coupling. 
Within this approach, only orientational spin fiuctua- 
tions are considered, whereas fluctuations of the ampli- 
tude of the local spin density are ignored. This approx- 
imation is justified below a crossover temperature Tx 
(of the order of the mean-field AF transition temper- 
ature) where the AF correlation length starts to grow 
exponentially (renormalized classical regime) . The effec- 
tive action of spin fiuctuations is then given by a NLctM. 
Solving the NLcrM within a "large-A/"" approach, we find 
that the ground-state of the Hubbard model on a square 
lattice is antiferromagnetic (Neel order) for any value of 
the Coulomb interaction U (Fig. 1).lI 

We have obtained the fermion spectral function 
yl(k, Lo) in the weak-coupling limit by computing the self- 
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energy E(k, w) to lowest order in the spin-fermion inter- 
action (Fig. 2). The QP peak which characterizes the 
Fermi Hquid state is suppressed by spin fluctuation when 
T <^ Tx- Instead, A(k, tj) exhibits a pseudogap separat- 
ing two broadened peaks. These peaks are precursors of 
the Bogoliubov QP's that appear at the T — Q AF tran- 
sition. Our results are in very scwd agreement with those 
obtained by the TPSC theory.El'aEl An important limita- 
tion of our analysis comes from the large-A/" solution of 
the NLctM. The latter introduces amplitude fluctuations 
of the Neel field which should be frozen at low tempera- 
ture. As a result, when going beyond the lowest-order 
contribution to S(k, w), we do not obtain the correct 
T limit of the fermion spectral function. In Ref. |l^, 
we show how this difficulty can be circumvented. 

There are several directions in which this work could 
be further developed. Since the NLcrM description is 
valid both at weak {U ^ t) and strong {U ^ t) cou- 
pling, our analysis of the fermion spectral function could 
be extended in the regime U ^ t. In the Mott-Hubbard 
insulator, we expect the pseudogap to transform into a 
(charge) gap of order U, the precursors of the T = AF 
bands becoming the upper and lower Hubbard bands. 

It is also possible to consider variants of the square lat- 
tice Hubbard model [Eq. where antiferromagnetism 
becomes frustrated. This would be the case for the t — t' 
Hubbard model (t' is the hopping amplitude for next- 
nearest neighbors) or if the lattice is triangular instead 
of square. Doping may also induce some kind of magnetic 
frustration.Ej This opens up the possibility to reach the 
quantum disordered and quantum critical regimes of the 
NLcrM (Fig. 1) and to study the corresponding fermion 
spectral functions. 
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found in the literature: incommensurate spin-density wave, 
spiral phase, stripes, etc. 
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FIG. 1. Crossover diagram derived from the large- A/" limit 
of the 2D NLcrM. At T = 0, there is long-range order when 
the dimensionless coupling g < Qc ~ 4:Tt. The three fi- 
/ ^ nite-temperature regimes correspond to "renormalized clas- 

/ sical" (RC), peiuantum critical" (QC) and "quantum disor- 

dered" (QD)y The ground-state of the half-filled 2D Hub- 
/ bard model on a square lattice is ordered for any value of the 

R-C ~~ ^ / QD Coulomb repulsion U . At finite temperature, there are strong 

\ ^ ^ AF fluctuations with an exponentially large coherence length 

^ / (RC regime) . 



QC 
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Hubbard model 
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FIG. 2. Lowest-order contribution to the fermion 
self-energy S. The dashed line represents the spin propagator 



X [Eq. (3.6) 
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